disp('1 vs 2');
model = fisher_2(rgb1, rgb2);
projC1 = linproj(rgb1', model);
projC2 = linproj(rgb2', model);
projCTest = linproj(rgb1_t', model);
[failure success] = bayes_2(projC1', projC2', projCTest', 1);

model = fisher_2(rgb1, rgb2);
projC1 = linproj(rgb1', model);
projC2 = linproj(rgb2', model);
projCTest = linproj(rgb2_t', model);
[failure success] = bayes_2(projC1', projC2', projCTest', 2);



str = sprintf('\n1 vs 3');
disp(str);
model = fisher_2(rgb1, rgb3);
projC1 = linproj(rgb1', model);
projC2 = linproj(rgb3', model);
projCTest = linproj(rgb1_t', model);
[failure success] = bayes_2(projC1', projC2', projCTest', 1);

model = fisher_2(rgb1, rgb3);
projC1 = linproj(rgb1', model);
projC2 = linproj(rgb3', model);
projCTest = linproj(rgb3_t', model);
[failure success] = bayes_2(projC1', projC2', projCTest', 2);



str = sprintf('\n1 vs 4');
disp(str);
model = fisher_2(rgb1, rgb4);
projC1 = linproj(rgb1', model);
projC2 = linproj(rgb4', model);
projCTest = linproj(rgb1_t', model);
[failure success] = bayes_2(projC1', projC2', projCTest', 1);

model = fisher_2(rgb1, rgb4);
projC1 = linproj(rgb1', model);
projC2 = linproj(rgb4', model);
projCTest = linproj(rgb4_t', model);
[failure success] = bayes_2(projC1', projC2', projCTest', 2);



str = sprintf('\n2 vs 3');
disp(str);
model = fisher_2(rgb2, rgb3);
projC1 = linproj(rgb2', model);
projC2 = linproj(rgb3', model);
projCTest = linproj(rgb2_t', model);
[failure success] = bayes_2(projC1', projC2', projCTest', 1);

model = fisher_2(rgb2, rgb3);
projC1 = linproj(rgb2', model);
projC2 = linproj(rgb3', model);
projCTest = linproj(rgb3_t', model);
[failure success] = bayes_2(projC1', projC2', projCTest', 2);



str = sprintf('\n2 vs 4');
disp(str);
model = fisher_2(rgb2, rgb4);
projC1 = linproj(rgb2', model);
projC2 = linproj(rgb4', model);
projCTest = linproj(rgb2_t', model);
[failure success] = bayes_2(projC1', projC2', projCTest', 1);

model = fisher_2(rgb2, rgb4);
projC1 = linproj(rgb2', model);
projC2 = linproj(rgb4', model);
projCTest = linproj(rgb4_t', model);
[failure success] = bayes_2(projC1', projC2', projCTest', 2);



str = sprintf('\n3 vs 4');
disp(str);
model = fisher_2(rgb3, rgb4);
projC1 = linproj(rgb3', model);
projC2 = linproj(rgb4', model);
projCTest = linproj(rgb3_t', model);
[failure success] = bayes_2(projC1', projC2', projCTest', 1);

model = fisher_2(rgb3, rgb4);
projC1 = linproj(rgb3', model);
projC2 = linproj(rgb4', model);
projCTest = linproj(rgb4_t', model);
[failure success] = bayes_2(projC1', projC2', projCTest', 2);